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Transformation-based cylindrical cloaks and concentrators are illuminated with non-monochromatic waves 
and unusual effects are observed with interesting potential applications. The transient responses of the devices 
are studied numerically with the Finite-Difference Time-Domain method and the results are verified with an- 
alytical formulas. We compute the effective bandwidth of several cloaking schemes as well as the effect of 
losses on the performance of the structures. We also find that narrowband behavior, frequency shift effects, time 
delays and spatial disturbances of the incoming waves are dominant due to the inherently dispersive nature of 
the devices. These effects are important and should be taken into account when designing metamaterial-based 
devices. 



I. INTRODUCTION 

Transformation electromagnetics provides the unprecedented opportunity to manipulate the propagation of electromagnetic 
waves in an artificial way [[1112] • Form-invariant coordinate transformations cause the electromagnetic waves to perceive space 
compressed or dilated in different coordinate directions, which often requires materials with anisotropic and spatially-varying 
values of permittivity and permeability tensors. Many interesting applications have been proposed, derived from different co- 
ordinate transformations. One of the most widely studied application is the cloak of invisibility [1], an approximate design of 
which was constructed for microwaves using arrays of split-ring resonators [3]. Other proposed devices include carpet cloaking 
structures lA^S], "cloaking" absorbers [6], concave mirrors for all angles [7], field concentrators [8], spherical [9] and cylindri- 
cal ifioll superlenses, flat near-field and far-field focusing lenses [11], reflectionless beam shifters and beam splitters lfl2h . field 
rotators ifTsI], adaptive beam benders and expanders 1 14] and novel antenna designs |[T5l[T6ll . 

In order to achieve the required material values, the majority of the aforementioned devices are constructed with metamaterials 
consisting of resonating structures [Jjl ll8|], whose properties are usually dominated by high loss and frequency dispersion 
im . Any material that exhibits permittivity and permeability values smaller than unity must be dispersive due to causality 
constraints ll2l|] . and hence, the corresponding devices operate properly only over very narrow bandwidths. One common 
example of such dispersive materials found in nature is the response of plasma electron gas to external radiation [22], which 
exhibits an effective refractive index smaller than one. Consequently, unless such a metamaterial-based device is illuminated by 
a unrealistic perfect monochromatic wave, a complete description of the response of such devices requires consideration of the 
full dispersive effects occurring over a range of frequencies. The goal of this paper is to investigate the physics of the broadband 
response of such dispersive devices. 

So far, two-dimensional (2D) cylindrical electromagnetic cloaks have been studied only under monochromatic plane wave il- 
lumination ill H 123,12^126]. However, all electromagnetic waves extend over a finite bandwidth. Some dispersive effects have 
been revealed analytically for three-dimensional (3D) spherical cloaks Jlzl Hsl HHl- For example, blueshift effects have been 
theoretically predicted |27], where waves with different frequencies penetrate differently inside the cloaking shell region. Fur- 
thermore, peculiar energy transport velocity distributions inside the same device have been predicted |28], using a Hamiltonian 
optics approach. Recently, the spatial energy distribution of a Gaussian light pulse after passing through such a spherical cloak 
has been shown to be distorted in a theoretical full- wave analysis (29] . However, these interesting effects have not been verified, 
until now, with a fully explicit numerical technique and especially for 2D cylindrical cloaks. In addition, the investigation of 
the cloaking bandwidth, which is constrained by the resonant nature of the metamaterial elements, has been very limited in the 
literature to mostly analytical treatments |27, 30]. 

In this paper, we utilize numerical simulations to examine the physics behind the dispersive and transient nature of 2D lossless 
cylindrical cloaks. The ideal |24] and matched reduced 1 3 1] cylindrical cloaks are illuminated with non-monochromatic Gaussian 
electromagnetic pulses, and their transient response is analyzed with a radially-dependent dispersive Finite-Difference Time- 
Domain (FDTD) method |26, 32]. The time-domain numerical technique presented here is advantageous compared to the Finite 
Element Method (FEM) used in previous works |24], since the transient response and the operational bandwidth of a device 
can be easily computed. First, we take advantage of the time-domain properties of FDTD to compute the bandwidth of some 
popular cloaking schemes. We show that the cloaking performance is very limited when reasonable losses are included in the 
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model. We also find that, under such non-monochromatic illumination, even the ideal electromagnetic cloak can be detected due 
to significant frequency shifts that occur because of the dispersion in the cloaking shell. This can be explained by the fact that the 
behavior of the waves at frequencies below the operating frequency of the cloak, where a negative-index shell may be formed, is 
fundamentally different than the behavior of waves above that frequency, where the cloak is found to act more like a dielectric 
scatterer. The Drude dispersion model is used to explain the origin of these results. In addition, variations in the group velocity 
and energy distribution across dispersive materials are found to cause severe distortion of the recomposed wavefronts after the 
cylindrical cloak and the electromagnetic concentrator devices. As a result, the incident pulses are distorted both temporally as 
well as spatially. These findings are not limited to the transformation-based cloaks examined in this work, but are expected to 
arise in a similar fashion in any dispersive and anisotropic device based on transformation electromagnetics. 

The rest of the paper is organized as follows. In Sec. Ullwe present the details of the dispersive FDTD code utilized in the 
numerical simulations. In Sec. [1111 the bandwidth and loss limitations of different cylindrical cloaks are explored. In Sec. |IVl 
the dispersive effects in the ideal cylindrical cloak are analyzed by launching non-monochromatic pulses towards the device, 
where frequency shift effects are examined. In Sec. [Vlthe corresponding dispersive effects in the reduced cylindrical cloak are 
investigated. Finally, in Sec. [Vllthe transient responses of the ideal cylindrical cloak and the ideal concentrator are investigated, 
which give rise to time-delay effects and non-uniform transmitted spatial energy distributions. 



II. RADIALLY-DEPENDENT DISPERSIVE FDTD 



In this section details of the FDTD simulations relevant to the modeling of dispersive devices are presented. For the 2D FDTD 
simulations, transverse electric (TE) polarized wave incidence is assumed, without loss of generality, where only three field 
components are non-zero: Ex^ Ey and Hz. The FDTD cell size, throughout the modeling, is chosen Ax = A?/ = A/ 150, where 
A is the wavelength of the excitation signal in free space. The domain size is 850 x 850 cells, or approximately 5. 66 A x 5. 66 A. The 
temporal discretization is chosen according to the Courant stability condition [33] and the time step is given by At = Ax/a/2c, 
where c is the speed of light in free space. Throughout the paper the devices are designed to have an operating frequency of 
/o = 2.0 GHz, where the free space wavelength is A :^ 15 cm. The cloaked object is chosen to be a perfect electric conductor 
(PEC) material. Unless otherwise noted, the dimensions of the cloaking structure are Ri = and = ^ in terms of the 
free space wavelength. Here Ri is the inner radius (radius of PEC cylinder) and R2 the outer radius of the cloaking shell. The 
FDTD computational domain used throughout this paper is shown in Fig. [Ha). 

Plane waves centered at different frequencies with variable temporal narrowband and broadband Gaussian envelopes are used 
to illuminate the structures. The Hz field amplitude values after the cloak are averaged along two overlapping line segments with 
different lengths, Li and L2, close to the right side of the domain, as can be seen in Fig. [Ha). These segments are located at the 
same point along the y-axis, despite being depicted at slightly different locations for clarity purposes. The frequency spectra of 
the transmitted pulses are then retrieved after recording the time history of the evolution of the field amplitude along these line 
segments. 

In this paper, the two main examples of interest are the ideal cloak and the matched reduced cloak. The parameters of the 
former device are given in cylindrical coordinates as 1I24I1 : 
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The parameters of the matched reduced cloak are instead 1I3TI1 : 
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Here Ri < r < R2is an arbitrary radius inside the cloaking shell. 

The material parameters of the devices discussed here can have values less than one, as extracted from the coordinate trans- 
formations, indicating that the required materials are frequency dispersive. The dispersive region of the parameters is mapped 
with the Drude model 1341 in the FDTD technique. For the results shown in this paper, it should be noted that physically similar 




FIG. 1 : (a) FDTD computational domain of the cylindrical cloaking structure for the case of non-monochromatic plane wave illumination. 
The fields are computed using the Total-Field Scattered-Field (TF-SF) technique | 33]. (b) Comparison between Drude and Lorentz dispersion 
models, when used to map the ir parameter at the point r = Ri. The bandwidth of a typical incident pulse is also plotted in the same graph 
to illustrate the dispersion along its spectrum. 



results are obtained, when the simulations are repeated using the Lorentz dispersion model (v?]. For the case of the radial 
components of the frequency-dependent relative permittivity ir and permeability fir in the FDTD code, they are given by: 

er[r,u;) = l ^ (3) 



uo'^ (r) 

where uops and uop^ are the plasma frequencies and 7 is the collision frequency, which is set to zero for the current lossless cloaks. 
The plasma frequencies are equal to: 



UJpeir) = - €r{r) (5) 



^Pfi{r) = iJoy^l - iir{r) (6) 

Here ujq is the device's operating frequency, while er{r) and /ir (^) are the desired frequency-independent material parameters 
at the design frequency cjq of the device. Equations (O and © are substituted in Eq. © and the resulted material parameters of 
the cloak are explicitly obtained: 

6r(r,u;) = 1-4[1- ^r(r)] (7) 



.(r,c.) = l-^[l-Ai,(r)] 



(8) 



Note that the model is build such that the device's design parameters are retrieved exactly when lo = luq. In that case Srir, toa) = 
er{r) and fir{r,u}o) = /Lir(r). 



Figure [TJb) shows the value of the real part of Sr as a function of frequency at the inner surface of an ideal cylindrical cloak, 
when either the Drude or Lorentz models are used. As a reference, the spectrum of an incident Gaussian pulse centered at 
2.0 GHz with a bandwidth of 200 MHz (Full Width at Half Maximum - FWHM) is also plotted in the same graph, in order 
to illustrate the variation due to frequency dispersion of the material parameters across a typical pulse. The graph shows that, 
independently of which dispersive model is used, waves at some frequencies perceive the inner core of the cloak as a negative- 
permittivity material, while waves at other frequency bands perceive it as a positive-permittivity material. As it shall be shown, 
this asymmetry in frequency, which is inherent to all dispersive materials, is the fundamental reason that gives rise to the effects 
that are described in the sections HVl and IVl 

The devices under consideration are anisotropic and highly dispersive. It has been shown in llssll that spatial resolutions of 
Ax < A/80 are necessary for the FDTD modeling of dispersive left-handed media, in order to avoid spurious resonances which 
are caused by numerical errors. Moreover, there are discrepancies between the analytical and the numerical values of the material 
parameters due to the finite spatial grid size used in FDTD simulations. In earlier work on the computational aspects of cloaking 
devices [26], we have demonstrated that spatial resolution smaller than A/80 is also necessary to accurately describe the material 
parameters in such schemes. Hence, in the work presented in this paper, a very fine spatial resolution of Ax = A/ 150 is safely 
chosen in order to avoid such issues. 

Another major challenge in FDTD simulations of dispersive devices based on transformation electromagnetics is late-time 
numerical instabilities [32]. These instabilities are mainly caused from the strong variations of the spatially varying material 
parameters. The extreme parameter values, combined with the staircase approximation of the device's cylindrical structure, 
can lead to spurious cavity resonances, which are visualized as accumulated charges at the interfaces between the boundaries 
of the device and the surrounding space. In order to mediate these instabilities, we have applied a locally spatial averaging 
technique | 36] for the anisotropic field components of the constitutive equations. Finally, corrected values of the plasma and 
collision frequencies have been used throughout the dispersive simulations, which are taking into account the finite time step of 
the FDTD technique. 

III. BANDWIDTH OF DISPERSIVE CYLINDRICAL CLOAKS 

In this section we use the time-domain capabilities of the FDTD technique in order to evaluate the bandwidth of certain 
cloaking devices. Since dispersive devices need to be constructed from resonant metamaterial elements in order to achieve the 
designed non-conventional material parameters, the cloaks theoretically operate correctly only for the single frequency for which 
the parameters of Eqs. ([B or Q are satisfied. Here we define the cloak's bandwidth through the frequency range over which 
the transmitted field components are enhanced compared to the corresponding transmitted field components observed when the 
same excitation input pulse is impinging on a bare non-cloaked PEC cylinder. 

The domain of the FDTD simulations can be seen in Fig. [Ha) and the dimensions of the devices are the same as mentioned in 
section ini The devices are excited with a plane wave pulse centered at 2 GHz confined in a broadband Gaussian envelope with 
a FWHM bandwidth of 1 GHz. The magnetic field values are spatially averaged along the parallel to the x-axis line segment 
1/2, approximately 2. 8 A away from the device's core. The transmitted spectrum is then retrieved from the time-dependence of 
the averaged field signals, which is next divided by the spectrum of the input pulse, yielding the transmission amplitude as a 
function of the frequency. 

The bandwidth performance of the ideal cylindrical cloak [24], the matched reduced cylindrical cloak [31] and the practical 
reduced cylindrical cloak [37] are compared to the transmission amplitude when a bare PEC cylinder is illuminated. The latter 
provides a reference point as we assume that a cloak operates only when it enhances the transmission amplitude compared to the 
bare cylinder case. The comparison between the transmission amplitudes of the different cloaking devices is shown in Fig. [2ja). 

We observe that the ideal cloak exhibits a transmission amplitude of 1 at the operating frequency of 2 GHz, when the corre- 
sponding value in the case of the bare cylinder is approximately 0.8. The calculated effective bandwidth is 11.5%. Among the 
proposed cloaks, this is the most broadband device as the bandwidths of the practical reduced cloak and the matched reduced 
cloak are only 4.6% and 2.5%, respectively. Thus, the approximate nature of the latter two devices significantly impacts its 
bandwidth performance. As expected, the transmission amplitude of both approximate devices is lower compared to the ampli- 
tude of the ideal cloak. However, the scattering performance of the matched reduced cloak is better compared to the practical 
reduced cloak, because its design allows it to be matched to the surrounding free space. Thus, there is a trade-off between the 
bandwidth and the scattering performance for the two approximate cloaking designs. 

As a final note, the cloaking bandwidth depends on the exact design of the cloak and its corresponding material parameters. 
For example, the PEC cylinder in these scenarios can be covered with thinner or thicker cloaks, other than the choices we have 
made here. Thinner cloaks require more extreme material parameters, as indicated by Eqs. ([B and ©, and are thus expected to 
be even more narrowband. On the other hand, by increasing the cloak's dimensions, its effective bandwidth is also expected to 
increase. 

One of the major practical issues with metamaterial-based devices is the effect of the losses that are inherent in the dispersive 
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FIG. 2: (a) The transmission amplitude of the ideal, matched reduced and practical reduced cylindrical cloak. The transmission of a bare 
PEC cylinder is also plotted to calculate the cloak's effective bandwidth, (b) The transmission amplitude of the ideal cylindrical cloak when 
different losses are introduced. 



resonating elements required to obtain the design material parameters ||2Q|] . which may inhibit their proper operation. Here 
we test the effect of losses in the cloaking bandwidth of the ideal cylindrical cloak. Similar as earlier in this section, FDTD 
simulations are performed with identical simulation domain and dimensions of the structure. Losses are introduced in the FDTD 
code in the same way as in our previous works |J, [12], i.e. by including a complex component — jtan(5 into the material 
parameters, which is increased up to 0.1. 

The transmission amplitude of the ideal cylindrical cloak as a function of frequency for various loss tangents can be seen 
in Fig. [2tb). We observe that the performance of the device is significantly affected with the increase of the losses, as higher 
losses decrease the transmission amplitude as well as the bandwidth of the cloak. It is interesting that if tan S > 0.05, which is 
a typical loss value for metamaterial structures close to resonance, then the performance of the ideal cloak is worst than a bare 
PEC cylinder and the cloaking phenomenon effectively ceases to exist. Similar results are observed when losses are introduced 
into the approximate cloaking structures (not shown), which implies that a practical, real world cloak will eventually require 
very low-loss metamaterial elements in order to operate properly. 



IV. SPECTRAL RESPONSE OF THE IDEAL CYLINDRICAL CLOAK 



In this section, the spectral response of the lossless ideal cylindrical cloak is investigated. The ideal cloak is excited with a 
narrowband Gaussian pulse with a bandwidth of 200 MHz (FWHM), centered at a frequency of /o = 2.0 GHz. The pulse is 
chosen to be narrowband in order to pass largely unaffected through the ideal cloak's allowed bandwidth (see Fig. [2la)). The 
results of the spectral content of the transmitted pulses recorded on each line segment of Fig. [Ha) are seen in Fig. [S] which are 
compared to the spectrum of the incident pulse. 

We observe that the central frequency of the transmitted pulse is shifted significantly only for the averaged field values 
corresponding to the shorter line segment Li, where the field distribution is strongly affected from the presence of the cylindrical 
cloak. When averaging over the longer segment L2 most of the pulse is transmitted unaffected and, as a result, the frequency 
shift effect is not observed. Only a small narrowing of the pulse is detected when averaging over L2, which is a result of the 
finite bandwidth of the device (Fig. [2ta)). For this simulation scenario, the blueshift is calculated to be A/ = 22.5 MHz, which 



FIG. 3: Blueshift effect observed in the normalized frequency spectra of transmitted Gaussian narrowband pulses through the ideal cylindrical 
cloak for two averaging line segments Li and L2 of Fig. \\Isl). The line segments are on the same position on the y-axis, 1.5A away from the 
cloak's outer shell, but have different lengths along the x-axis (Li — Ri, L2 — 4.25i?2)- The effect is stronger for the line segment Li, near 
the center of the cloaking structure. 

is 1.1% deviation from the central frequency of 2 GHz. Hence, an instrument which is capable of resolving spectral deviations 
smaller than 1.1%, will detect the presence of the cloak. 

To further investigate the blueshift effect, six narrowband Gaussian pulses are independently launched towards the ideal 2 
GHz cloak, at central frequencies of 1.7, 1.8, 1.9, 2.1, 2.2, 2.3 GHz respectively. The FWHM bandwidth of each pulse is 200 
MHz. The spectral distribution of the averaged field values of the transmitted pulses along the line segments Li and L2 depicted 
in Fig. (Ha) 1.5 A away of the cloak is once again evaluated, and is shown in Figs. Oc), (d). In addition, we also record the 
transmitted spectral distribution along two new line segments of the same length Li and L2 using the same excitation signals, 
but this time the segments are positioned immediately after the cloaking structure. These results are shown in Figs. (Ha), (b). 

Firstly, it is observed that when the fields are averaged over the L2 fine (Figs. Htb), (d)) the differences between the pulses 
above and below the cloak's central frequency are small. The spectrum is affected symmetrically around the central frequency 
because of contributions from field energy away from the device that does not interact with the cloak. The lower amplitudes 
recorded away from the central frequency region are caused by the overall finite bandwidth of the cloak, as shown in Fig. [SJa). 
Secondly, when recording over the short line segment Li (Figs. HJa), (c)), it is observed that the pulses with central frequencies 
greater than the cloak's operating frequency (/ > /o) maintain their overall shape but are reinforced in amplitude. On the 
contrary, the pulses with central frequencies less than the cloak's operating frequency (/ < /o) are dissipated, when they are 
passing through the structure. This behavior gives rise to the frequency shift effect observed in Fig. [3l The effect is stronger 
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FIG. 4: Frequency spectra of six transmitted narrowband pulses, centered at different frequencies, after impinging on the ideal cloak. The 
cloak's design frequency is 2 GHz. The fields are averaged over two different line segments Li ((a), (c)) and L2 ((b), (d)) which are perpen- 
dicular to the direction of propagation as shown in Fig. [ita). The line segments are positioned either right after the cloak's outer boundary 
((a), (b)) or at a distance 1.5 A away ((c), (d)). 
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FIG. 5: (a) Analytical normalized penetration depths for ideal cylindrical cloaks with different dimensions, (b) Regions formed inside the 
ideal cylindrical cloak, when it operates at / < /o frequencies, (c) Analytical normalized penetration depths for matched reduced cylindrical 
cloaks with different dimensions. The numerically estimated penetration depths are also shown for a cloak with thickness 2 A/3. The error bars 
indicate numerical uncertainty on the field cutoff radius, (d) Regions formed inside the matched reduced cylindrical cloak, when it operates at 
f < fo frequencies. 



away from the cloak (Fig. IHc)), as opposed to near the cloak (Fig. lUa)), since the field has recombined behind the device 
in the former case. At an even greater distance away from the device, the amplitudes of the transmitted pulses do not evolve 
significantly. 

We now investigate the effective structure that an off-frequency incident wave perceives inside the dispersive cloaking shell. 
Depending on the frequency of the incident wave, the material parameters of the ideal cloak (Eq. ([T])) may become equal to 
zero somewhere inside the cloaking shell. These locations are then effectively acting as a PEC wall (when £r = 0) and perfect 
magnetic conductor (PMC) wall (when /i^ = 0) for that frequency, beyond which the fields theoretically do not penetrate |27]. 
The radial locations r^, r^, where the permittivity and permeability respectively vanish for the ideal cloak, are given as a function 
of frequency by (using Eqs. (|7]), ([5]), O): 

^2 

r^(c^<c^o) = — -2 

The theoretical field penetration depths normalized to the object's radius R\ are plotted in Fig. Ha) for different dimensions 
of the cloak as a function of frequency, for f < fo (for / > /o no such walls are formed). It is observed that always > r^. 
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FIG. 6: Frequency spectra of six transmitted narrowband pulses, centered at different frequencies, after impinging on the matched reduced 
cloak. The cloak's design frequency is 2 GHz. The fields are averaged over two different line segments Li (a) and L2 (b), as shown in Fig. 
[U^a), which are perpendicular to the direction of propagation. The line segments are positioned at a distance 1.5A away from the cloak's outer 
shell. 

with Vs for any cc; as ^ ^1 • This implies that a PEC wall will form first as the wave approaches the core of the 

device, beyond which < 0. If a material with positive index of refraction was formed between that wall and the device's 
PEC core, no fields would penetrate further inside. However, further inside the cloak a PMC wall will be then formed, beyond 
which also /i^ < 0. Thus, a negative index shell will be formed for some frequencies near the core of the device, which can 
trigger phenomena such as superscattering |7]. This three-layer effect is unique to the ideal cylindrical cloak (for TE or TM 
polarizations) because two material parameters can become simultaneously negative due to dispersion, something that is not 
occurring in either the ideal spherical cloak or the matched reduced cylindrical cloak. These layers are depicted graphically in 
Fig. Ob) for / < /o frequencies. 

In addition, since for f = fo the material parameters become zero at the inner core r = Ri, waves with frequencies f < fo 
will sample negative values of the material parameters near the core, while waves with frequencies f > fo will sample strictly 
positive material parameters values everywhere in the cloaking shell. Thus, the former waves (/ < /o) will perceive the device 
as conducting scatterer, while the latter (/ > fo) will perceive it as a dielectric and magnetic material, with positive permittivity 
and permeability. 

This latter effect, that waves at different frequencies perceive electromagnetically different devices, is demonstrated by launch- 
ing three monochromatic plane waves at frequencies 1.7 GHz, 2.0 GHz and 2.3 GHz against the ideal cloak. The real part of 
the converged magnetic field values are shown in Fig. |3a)-(c). We observe in Fig. |7ta) that for the 1.7 GHz wave, little field 
reaches the core of the device due to the aforementioned walls that form, and thus a strong shadow is observed behind the device. 
On the contrary, in Fig. |3c), the 2.3 GHz wave perceives the cloaking shell as a positive index material and field is distributed 
even inside the device. As a reference, the field distribution on the nominal frequency of 2.0 GHz is also shown in Fig. |7lb). 
The positive-to-negative material parameters transitions that occur for f < fo but not for f > fo give rise to the frequency 
modulation that occurs inside the cloaking shell (shown in Fig. lUa), (c)), as waves at / < /o are dissipated due to the scattering 
from a conductor, while others at / > /o are scattered from a positive-index material. This ultimately produces the frequency 
shift observed in Fig. [3] 



We now focus on the matched reduced cylindrical cloak L31J . which is easier to be practically implemented compared to the 
ideal cloak. This material parameter set is more similar to the ideal spherical cloak JTI], in the sense that, for a given polarization, 
only one parameter is radially-dependent and dispersive (£r for TE waves), whereas the others are constant and conventional. 
FDTD simulations are performed over the same domain as before (Fig. [Ha)) and the cloak has the same dimensions. Again, 
narrowband Gaussian pulses at various central frequencies are launched towards the cloak. The transmitted pulses in the fre- 
quency domain, after impinging on the cloak, recorded at the two line segments Li and L2 of Fig. [Ha), at a distance 1.5 A away 
from the device, are seen in Figs. Oa) and (b), respectively. 

The transmitted pulses appear slightly distorted compared to the ideal cloak (Fig. [4]) due to the fact that the reduced cloak ex- 
amined here is designed to operate approximately. However, the features of the transmitted pulses are similar to the transmission 
observed through the ideal cloak, where frequencies f < fo are dissipated while frequencies f > fo are enhanced. The effect is 
once again much stronger when averaging the fields near the cloak's core, and gives rise to a similar frequency blueshift effect. 

For some frequencies / < /o, the dispersive parameter ir will become equal to zero at a specific location inside the cloaking 
shell, forming a boundary beyond which wave penetration is very weak. The penetration depth is calculated analytically with a 



V. SPECTRAL RESPONSE OF THE REDUCED CYLINDRICAL CLOAK 
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FIG. 7: Real part of the magnetic field amplitude distributions (Hz) when plane waves of different frequencies are impinging on the ideal 
((a)-(c)) and matched reduced ((d)-(f)) cloak, after steady state is reached. The amplitude scale is normalized such that 1 is the maximum 
plane wave field amplitude without any device present. (a),(d): At 1.7 GHz a PEC wall is formed inside the cloaking shell where fields cannot 
penetrate. The devices behave as a conductive scatterer and a large shadow is observed behind them. (b),(e): At 2.0 GHz the devices operate at 
their nominal frequency. For (e) only, imperfections in the field distributions are inherent in the limitations of the cloak's approximate design. 
(c),(f): At 2.3 GHz the cloaking material is perceived as a dielectric scatterer from the incident wave. Thus no significant shadow is formed 
behind the devices. 



similar procedure as before and it is given by the equation: 

r,{u; < coo) = , — (11) 

It is normalized to Ri and plotted as a function of the / < /o frequencies in Fig. Oc) for cloaks with different thicknesses. 
Moreover, the locations, where the impenetrable wall is formed, are deduced approximately through the FDTD simulation by 
identifying the points where the magnetic field Hz becomes zero. They are also shown in Fig. Oc) for a specific cloak design 
and are found to be in excellent agreement with the analytical predictions of the penetration depth derived from Eq. (fTTI) . The 
layers that form inside the reduced cloak are depicted graphically in Fig. [5td) for / < /o frequencies. 

Similar to the procedure followed in the previous section for the ideal cloak, in order to explain the different behavior of the 
pulses above and below the central frequency, along with the formation of the PEC wall, three monochromatic plane waves are 
launched with frequencies 1.7, 2.0 and 2.3 GHz towards the 2 GHz matched reduced cloak. The real part of the converged field 
amplitude values are seen in Figs. |7td)-(f). As a reference. Fig. (Tte) shows the field distribution at the cloak's nominal frequency, 
where the wave energy mostly recomposes smoothly behind the device. The imperfections observed in the wavefronts are due 
to the approximate nature of the reduced cloak's design. 

For the plane wave with frequency 1.7 GHz, according to Eq. (fTTI) . an impenetrable wall is expected to form approximately 
0.36A away from the inner core. This is indeed observed in Fig. Wid), where little field reaches the metallic core. As was 
discussed previously, due to the PEC wall the wave at that frequency should perceive the device as a conductive scatterer. This 
is confirmed in the simulation by observing the large shadow that appears behind the device, similar to the shadow that typically 
appears by the scattering of a plane wave off a conductive cylinder. This shadow is the reason for the reduced amplitude recorded 
for wave pulses at frequencies / < /o, as was pointed out in Fig. Oa). The shadow is similar to the one observed in the case of 
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FIG. 8: Time-dependent snapshots of the real part of the magnetic field amplitude distribution (Hz) when a 1 GHz (FWHM) Gaussian pulse 
is impinging on dispersive devices, (a)-(c): Ideal cloak, (d)-(f): Ideal concentrator. In (a),(d) the pulse has reached the first half of the devices. 
In (b),(e) the pulse leaves the dispersive regions. The wavefronts near the center of the devices are delayed compared to the wavefronts away 
from the central regions. In (c),(f) a delay appears near the center of the devices, even though the waves have recomposed. Reflections are 
observed due to the broad bandwidth of the incident pulse. 

the ideal cloak, as shown in Fig. [Ha). 

For the plane wave with frequency 2.3 GHz, no PEC wall is predicted to exist. Instead, the wave perceives the cloaking shell 
as a dielectric scatterer, since the relative permittivity values are always larger than zero for / > /o- Indeed, the field distribution 
shown in Fig. |3f) indicates that the fields penetrate all the way into the inner core, and that no shadow appears behind the 
device. On the contrary, the interference effects cause the field amplitude behind the cloak to be locally larger than the amplitude 
of the incident field (equal to unity in Fig. [7]), thus enhancing these / > /o frequencies compared to the / < /o frequencies. 
This explains the modulation pattern that appears in Figs. Oa) which ultimately causes the blueshift effect presented in Fig. [3] 

VI. TEMPORAL AND SPATIAL RESPONSES OF THE IDEAL CYLINDRICAL CLOAK AND CONCENTRATOR 

Another interesting effect takes place in the ideal cloaks for non-monochromatic waves. While light rays that travel away 
from the core of the cloak propagate at the speed of light in free space, rays that traverse the core region experience time delays 
[28], due to the increased values of the permittivities and permeabilities introduced by the anisotropic cloaking materials. This 
effect is stronger near the center of the cloak, where in theory some material parameters reach infinite (or realistically very large) 
values. As we shall show, the phenomenon is visible when a temporally finite pulse is incident towards the cloak: wavefronts 
propagating away from the core reach the other side of the domain (behind the cloak) sooner than wavefronts propagating near 
the center. 

The time-delay effect has been predicted theoretically for a 3D ideal spherical cloak Jlsl |29] . In this paper, this dynamic 
effect is demonstrated and exploited with the FDTD method for the ideal cylindrical cloak. The same dimensions, as before, are 
chosen for the tested cloak. A broadband Gaussian pulse is launched and is passing through the 2D structure. It has a bandwidth 
of 1 GHz, centered at the cloak's nominal frequency of 2 GHz. The real part of the magnetic field amplitude distribution is 
calculated as a function of time through the FDTD simulation. Three time snapshots of the pulse are depicted in Figs, [lla)-(c). 
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FIG. 9: Relative accumulated energy distribution for ideal cloak and concentrator. In the amplitude scale shown, 1 is the accumulated energy 
distribution of free space propagation. The lines Ri and R2 indicate the inner and outer radii of the devices. The distributions are not perfectly 
symmetric because of the finite discretization used in the FDTD simulations. 

as the pulse is propagating through the device. When the pulse is recomposed after the cloak (Fig. [Sc)), the wavefronts are not 
flat anymore due to the experienced time delay, which is more intense closer to the cloak's inner boundary. The group velocity 
is reduced close to the inner boundary of the structure, whereas the phase velocity is approaching large values. 

In this example, the wavefronts near the center of the cloak are delayed by approximately 0.27 ns compared to wavefronts 
that propagate mostly undisturbed away from the cloak. The bending of the waves is caused by the extremely high values of 
the radially-dependent cloaking parameters (Eq. ([T])) due to the dilation of the domain from a point to a circle ring, the device's 
inner core. Hence, the time delay of the pulse is inherent to the design of the cloak and, as a result, unavoidable. Thus, one could 
take advantage of this effect in order to detect the presence of the cloaking device with an appropriate instrument. It is present 
in both the ideal and reduced (not shown here) 2D cylindrical cloaks. It is also the main reason for the long time required for the 
FDTD simulations to reach steady-state results | 38] in such simulations. 

The time-delay effect is not limited to cloaks, but generally appears whenever large variations of the values of the material 
parameters are required for a given device. To illustrate a different scenario. Figs. [8td)-(f) show time snapshots as the same 
broadband Gaussian pulse is impinging on the ideal concentrator [8], a device designed to focus electromagnetic fields inside 
a small region of space (r < Ri in this case). The dimensions of the device are the same as the cloak's discussed earlier. The 
time-delay effect is observed after the wave recomposes behind the device, especially near the core, where the wavefronts are 
delayed by roughly 0.18 ns compared to the unperturbed waves. The effect is weaker compared to the cloak, since none of the 
material parameters reach extreme values (maximum value of permittivities and permeability is 3) | 8]. Meanwhile, the phase 
and group velocities have comparable values. It should also be noted that reflected waves emerge as the pulses impinge on either 
device (Figs. [He) and (f)). This is a result of the large number of frequency components contained in the Gaussian pulses, since 
the devices can only operate properly at their nominal frequency of 2 GHz. 

In addition to the temporal disturbances induced on the Gaussian pulses by the dispersive devices, their spatial distribution is 
affected as well ||2^ . To illustrate this effect, the time-integrated energy distribution crossing the line segment L2 (Fig. [Ha)) for 
both ideal cloak and concentrator is recorded in the setup of Fig. [8l and it is shown in Fig. [9j We observe that the energy of 
the non-monochromatic source is not distributed uniformly outside of the devices due to scattering of the frequency components 
away from /q. The cloak has more extreme parameters than the concentrator, which leads to less accumulated energy inside the 
device and more scattered energy outside. It is interesting that the energy reaching the center of both structures is lower than the 
scattered energy near the edges of the simulation domain. Hence, the energy is not properly spatially recombined behind either 
device. 



VII. CONCLUSIONS 

Concluding, the properties of the metamaterial dispersive devices under non-monochromatic illumination are studied. The 
inherent dispersive nature of the devices is the main reason behind the transient effects observed, which is affecting their broad- 
band performance. We investigate the origin of the frequency shifts that appear in both the ideal and matched reduced cylindrical 
cloaks through dispersive FDTD simulations, using narrowband Gaussian pulses as illumination. In addition, the time-delay and 
spatial non-uniformity effects that occur near the center of the dispersive regions of the ideal cloak as well as the ideal con- 
centrator are demonstrated. These are important effects that should be taken into account in future applications of dispersive 



metamaterial devices. 

The ideal cylindrical cloak and concentrator, work properly only for monochromatic wave incidence. However, they become 
non practical when they are excited with non-monochromatic radiation, which is the type of electromagnetic radiation found in 
nature. The effect of losses is also significant for such dispersive devices. In general, the transmitted signal, passing through 
devices derived from coordinate transformations, will have its spectral, temporal and spatial profiles altered. These restraints 
could be avoided if the cloak is constructed from broadband active metamaterials lUilli^]. Alternatively, novel non-dispersive 
devices based on conventional materials (47] need to be implemented. 
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